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Abstract. This paper predicts self-consistent faint 
galaxy counts from the UV to the submm wavelength 
range. The STARDUST spectral energy distributions de- 
scribed in Devriendt et al. ( 1999| ) (Paper I) are embedded 
within the explicit cosmological framework of a simple 
semi-analytic model of galaxy formation and evolution. 
We begin with a description of the non-dissipative and dis- 
sipative collapses of primordial perturbations, and plug in 
standard recipes for star formation, stellar evolution and 
feedback. We also model the absorption of starlight by 
dust and its re-processing in the IR and submm. We then 
build a class of models which capture the luminosity bud- 
get of the universe through faint galaxy counts and red- 
shift distributions in the whole wavelength range spanned 
by our spectra. In contrast with a rather stable behaviour 
in the optical and even in the far-IR, the submm counts 
are dramatically sensitive to variations in the cosmologi- 
cal parameters and changes in the star formation history. 
Faint submm counts are more easily accommodated within 
an open universe with a low value of Oo, or a flat universe 
with a non-zero cosmological constant. We confirm the 
suggestion of Guiderdoni et al. ( 1998 ) that matching the 
current multi-wavelength data requires a population of 
heavily-extinguished, massive galaxies with large star for- 
mation rates (~ 500 M© yr _1 ) at intermediate and high 
redshift (z > 1.5). Such a population of objects probably 
is the consequence of an increase of interaction and merg- 
ing activity at high redshift, but a realistic quantitative 
description can only be obtained through more detailed 
modelling of such processes. This study illustrates the 
implementation of multi-wavelength spectra into a semi- 
analytic model. In spite of its simplicity, it already pro- 
vides fair fits of the current data of faint counts, and a 
physically motivated way of interpolating and extrapolat- 
ing these data to other wavelengths and fainter flux levels. 

Key words: cosmology: structure formation - cosmology: 
galaxy counts 



1. Introduction 



In the dense gas clouds that harbour starbursts, the ultra- 
violet (UV) light of young stars is absorbed by dust grains 
which, in turn, release their thermal energy at infrared 
(IR) and submillimetre (submm) wavelengths. Thus, un- 
derstanding the star formation history of galaxies clearly 
requires a correct assessment of the UV to submm lumi- 
nosity budget. The most straightforward and simple obser- 
vational probe of such a luminosity budget is the analysis 
of the faint galaxy counts obtained at various wavelengths. 
In this purview, this paper proposes self-consistent theo- 
retical predictions of faint galaxy counts at optical, IR and 
submm wavelengths that can be directly compared with 
the current host of data, and used to prepare observational 
strategies with forthcoming instruments. 

In the local universe, only 30 % of the bolometric lu- 
minosity is released in the IR/submm wavelength range 
(Soifcr and Neugebauer, 1991), and the effect of extinc- 



tion can thus be considered as a mere correction that does 
not change the main evolutionary trends elaborated from 
optical studies. However, there is now a growing amount 
of evidence that this fraction was much higher in the past. 
Indeed, the discovery of the Cosmic Infrared Background 
(CIRB) at a level ten times higher than the no-evolution 
predictions based on the IRAS local IR luminosity func- 
tion, and twice as high as the Cosmic Optical Background 
obtained from optical counts, showed that dust extinc- 
tion and emission play a major role in defining the lu- 
minosity budget of high-redshift galaxies (Puget et al. 
( |l996| ), Guide rdoni et al. ( |l997| ), Schleg el et al. ( p^98| ), 
Fixsen et al. (1998), Hauser et al. ( 1998| )). Since this ma- 



jor breakthrough, deep surveys with the ISO satellite at 
15 /mi (Oliv er et al. ( |l997j ), Aussel et al. (|1999| ), Elbaz 
et al. ( |1999|) ) and 175 //m (Kawara et al. ( |1998 ), Puget 
et al. (199S)), and with the SCUBA inst rume nt at 850 
fim (Smail et al. ( |1997| ), Barger et al. ( |1998p , Hughe s 
et al. ( |1998| ), Eales et al. ( p^99| ), Barger et al. ( |l999a| )) 
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have begun to break the CIRB into its brightest contribu- 
tors. Although identification and spectroscopic follow-up 
of the submm sources are not easy, the preliminary results 
of such studies seem to show that part of these sources 
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are the high-redshift counterparts of the local luminous 
and ultraluminous IR galaxies (LIRGs and ULIRGs) dis- 
covered by I RAS (Sm ail et al. ( |1998| ), Lilly et al. ( |1999| ), 
Barger et al. (1999b)). In parallel to this pioneering explo- 
ration of the "optically-dark" and "infrared-bright" side 
of the universe, a more careful examination of the Canada- 
France Redshift Survey (CFRS) galaxies at z ~ 1, and 
Lyman Break Galaxies at z ~ 3 and 4 do show a signifi- 



ever, predictions at IR/submm wavelengths from a SAM 
were produced much later by Guiderdoni et al. ((1997), 



cant amo unt of extinction (Flo res et al. (1999), Steidel et 
al. (|1999| ), Meurer et al. (|l999|) ). The previous estimates 
of the UV fluxes, and consequently of the star formation 
rates, in these objects have to be respectively multiplied 
by factors 3 and 5 to take into account the effect of extinc- 
tion. Dust seems to be present at still higher redshif ts. For 
instan ce, it is seen in a lensed galaxy at z = 4.92 ( [goifcr 
et al., 1998) and even in a Lyman a galaxy at z — 6.68 
( |Chen et al., 1999| ). 

The synthetic spectra of stellar populations in galax- 
ies are easily computed from spectrophotometric models 
of galaxy evolution. Unfortunately, most of these models 
neglect the influence of dust on the spectral appea rance 
of galaxies. Guiderdoni & Rocca-Volmerange (1987) pro- 
posed a first modelling of the effect of dust extinction. 
Later, Mazzei et al. (1992) basically used the same recipe 
for extinction, but they also computed dust emission to 
get spectral energy distributions (SEDs) from the UV to 
the far-IR. Complete sets of synthetic spectra are now 
available from the GRASIL ( Silva et al., 1998| ) and star- 
dust models (Devriendt et al. (1999), hereafter Paper I) 
of spectrophotometric evolution. These models share the 
same spirit, but they differ by a number of details. The 
GRASIL SEDs in the IR are computed from a more sophis- 
ticated model of transfer, that is more explicit physically, 
but involves several free parameters, whereas Stardust 
SEDs in the IR are computed with a minimal number of 
free parameters, by weighing various dust components to 
reproduce the observed IRAS colour-luminosity relations. 

These spectra can be used in phenomenological mod- 
els of faint galaxy counts, that extrapolate the evolution of 
the local galaxies backwards under the assumption of pure 
luminosity evolution. For instance, the predictions of faint 
galaxy counts at o ptical wavelengths by Guiderdoni & 
Rocca-Volmerange (1990) used their optical spectra with 
extinction, whereas Franceschini et al. (1991), (1994) used 
the Mazzei et al. spectra with dust extinction and emission 
to produce the first set of counts at optical and FIR wave- 
lengths. However, semi-analytic models of galaxy forma- 
tion (hereafter SAMs) are a much more powerful approach 
to describe the physical processes that rule galaxy forma- 
tion and evolution within an explicit cosmological con- 
text (White fc Frc nk ( [l991~l ), La cey fc Silk ( |i99l| ), Kauff- 
mann et al. Q1993D, Cole et al. Ql994|), Somerville & Pri- 



mack ( |2000D ). Wh ite fc Frcnk ( |1991| ), Lac ey et al. ( |1993| ) 



Kauffmann et al. (1994), and Cole et al. (1994) proposed 



predictions of faint galaxy counts at optical wavelengths 
(basically the B and K bands) from their models. How- 



( |1998| ), hereafter GHBM). But this first study did not give 
the corresponding predictions at optical wavelengths, and 
was restricted to the = 1 standard Cold Dark Matter 
(CDM) model. 

In this paper, we implement the STARDUST spectra 
into a SAM to make predictions of faint galaxy counts 
and redshift distributions from the UV to the submm 
wavelength range, very much in the spirit of GHBM. We 
also extend the SAM to other cosmologies, and study the 
sensitivity of the results to the cosmological parameters 
and star formation history. Although our approach has a 
number of shortcomings which are due to the simplicity 
of our model, this paper primarily intends to show that 
(i) the implementation of STARDUST SEDs into SAMs is 
straightforward because of its small number of free pa- 
rameters; (ii) the model gives fits that are already very 
satisfactory in spite of the simplicity of the approach; and 
(iii) the outputs are a physically motivated tool to interpo- 
late or extrapolate the current observations of faint counts 
to other wavelengths and/or flux levels. This is needed to 
prepare the observational strategies with the forthcoming 
IR/submm satellites sirtf, first and planck, as well as 
with the Atacama Large Millimetre Array. 

In section 2, we briefly describe how we connect the 
Stardust spectra with the various physical processes that 
are relevant to galaxy formation, within our SAM. We 
point out the differences with GHBM. Section 3 discusses 
the values of the free parameters that define the so-called 
"quiescent" mode of star formation. Section 4 focuses on 
the sensitivity of the faint galaxy counts to a change in 
the cosmological parameters. Section 5 studies the sensi- 
tivity of the faint galaxy counts to galaxy evolution, and 
more specifically to the presence of a heavily-extinguished 
"starburst" mode of star formation similar to the one in 
ULIRGs. Finally, a fiducial model is proposed. We discuss 
our results in Section 6. 



2. The basics of our semi analytic model 

In the SAM ab initio approach, galaxies form from Gaus- 
sian random density fluctuations in the primordial matter 
distribution, dominated by CDM. Bound perturbations 
grow along with the expanding universe, until gravitation 
makes them turn around and (non-dissipatively) collapse. 
As a result, they end up as virialized halos. Then the 
collisionally-shocked baryonic gas cools down radiatively, 
and settles at the bottoms of the potential wells where it is 
rotationally-supported. Stars form from the cold gas and 
evolve. At the end of their lifetimes, they inject energy, gas 
and heavy elements back into the interstellar medium. The 
chemical evolution is computed, and the recipes developed 
in GHBM and in Paper I give the amount of optical lu- 
minosity that is absorbed by dust and thermally released 
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at IR/submm wavelengths. Finally, overall SEDs from the 
UV to the submm are computed with STARDUST. 

We refer the reader to GHBM for a detailed descrip- 
tion of how to compute the mass distribution of collapsed 
dark matter halos from the peaks formalism introduced 
by Bardeen et al. (1986), and Lacey and Silk (1991) in 



an Einstein-de Sitter universe. We give in appendix A the 
quantities which enable us to extend this formalism to low 
matter-density universes with or without a cosmological 
constant. As we follow closely the prescriptions in GHBM, 
we only mention in the following subsections the quanti- 
ties which differ from their work. 

2.1. Gas 

We assume that a universal "baryonic fraction" Os/fio of 
the pristine gas gets locked up within each dark matter 
halo, where it is collisionally ionised by the shocks occur- 
ring during virialization. Because it can cool radiatively, 
gas then sinks into the potential wells of the halos. The 
cooling time depends on the gas metallicity. Here, we de- 
cide to adop t the cooling function given by Sutherland 
and Dopita ( 1993 ) for one third of solar metallicity. This 
choice is motivated by the fact that it is the average value 
that is observed today in clusters, and probably, as argued 



by Renzini (1999), the average value of the low-redshift 
universe as a whole. Under this assumption, the cooling 
time is underestimated in high-redshift halos where the 
gas is more metal-poor. However, these objects are also 
smaller and denser on an average, so that their cooling 
times are already very short. 

We assume that the gas stops falling into the dark 
matter potential wells when it reaches rotational equilib- 
rium, and forms rotating thin disks (see e.g. Dalcanton et 
al. ( [1997D and Mo et al. fll998| )). Following these authors, 
we adopt for the thin disk an exponential surface density 
profile with scale length r^, and truncation radius r t , such 
as: 



E(r) = 



£(0)exp(-£) ifr<r* , 







if r > r t 



where r t is defined as the minimum value between the 
virial radius rw, and f c Td- The free parameter f c defines 
the extent of the cold gas disk. 

We then relate the exponential scale length of the cold 
gas disk ra. to the initial radius r V j r , through conservatio n 
of specific angular mome ntum ( Fall and Efstathiou, 1980| ) . 
As shown by Mo et al. ( 1998 ), stability criteria yield: 



f—s\t vtr 

y/2 



fc exp(-/ c ) 



2 [1 - (1 + f c ) exp(-/ c )] 



(1) 



where A = J\E\ 1 / 2 G- 1 M~ 5 / 2 ~ 0.05 ± 0.03 is the well- 
know dimensionless spin parameter. 

As this will be important later, we emphasise that 
the simple formalism used here does not allow us to form 



spheroids through mergers/interactions of galaxies. There- 
fore, in section 5, we will define a "starburst mode" which 
phenomenologically accounts for this process. 

2.2. Stars 

The only time scale available in our gas disks is the dy- 
namical time scale td yn = ^irr d/Vc- Therefore, guided by 
observational data (Kennicutt ( 1998) ) ) , we assume that the 
complicated physical processes ruling star formation lead, 
at least in a disk galaxy, to a global star formation rate 
(SFR) with the simple law: 



SFR(t) 



M gas (t) 



(2) 



where M gas (t) is the total mass of cold gas in the disk 
at time t. We introduce an efficiency factor /3 _1 as a sec- 
ond free parameter. The IMF is chosen to be Salpeter's, 
with slope x = 1.35 between masses md — 0.1 and 
m u = 120M Q . 

The Stardust spectrophotometric and chemical evo- 
lution model presented in Paper I is then used to compute 
metal enrichment of the gas as well as the UV to NIR 
spectra of the stellar populations produced with such star 
formation rates. Details on the stellar spectra, evolution- 
ary tracks and yields can be retrieved from this paper and 
references therein. 

2.3. Feedback 

Along with producing metals, massive stars which, at the 
end of their lifetimes, explode in galaxies, eject hot gas and 
heavy elements into the interstellar and/or intergalactic 
medium. We focus here on the modelling of this " stella r 
feedback", which is inspired from Dekel and Silk ( 1986 ). 
The average binding energy of a mass of gas M gas (t) 
distributed within a truncated exponential disk at time 
t, which is gravitationally-dominated by its dark matter 
halo, is given by: 

\ ^ 2nrM gas (r,tMr)dr ~ ^M gas (t) [V 2 + V 2 sc (r t )} 

where is the gravitational potential of a singular 

isothermal sphere truncated at virial radius 7W, and: 



V esc (r) = V c 



2 1 - In ■ 



"I 1/2 



(4) 



is its escape velocity at radius r. 

As a result, the energy balance between the gravita- 
tional binding energy and the kinetic energy pumped by 
supernovae into the interstellar medium yields the frac- 
tion of stars F± that formed before the triggering of the 
galactic wind (at time tw)' 



1 



M+(t w ) + M gas {t w ) 1 + {V hot /V c ) 2 ' 



(5) 
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with: 



Vhot = 



^VsnEsncsn 



1 + (V esc (r t )/V c y 



1/2 



(G) 



where the energy available per supernova is E$n = 
10 51 erg, and the number of supernovae per mass unit 
of the stars that just formed is tjsn = 7.5 x 10 -3 for 
a Salpeter IMF. The SN heating efficiency e$N is a third 
free parameter. By taking the initial gas mass available for 
star formation to be the initial cold gas mass minus the 
gas mass lost in the galactic wind, one then approximate 
chemical evolution by the closed-box model described in 
Paper I. Note that we have neglected any dynamical effect 
due to mass loss in the previous analysis. 

24. Dust 

Part of the luminosity released by stars is absorbed by dust 
and re-emitted in the IR/submm range. We now briefly 
outline how we compute the luminosity budget of our ob- 
jects within the SAM. We emphasise that this is a major 
improvement with respect to GHBM, as for the first time, 
stellar and dust emission are linked self-consistently. As 
in Paper I, we proceed to derive the IR/submm dust spec- 
tra with three steps: (i) computation of the optical depth 
of the disks, (ii) computation of the amount of bolomet- 
ric energy absorbed by dust, and (iii) computation of the 
spectral energy distribution of dust emission. 

The first step is easily completed because we know 
the sizes of our objects from eq. [l] and the definition 
of the truncation radius r t , and we obtain the mass of 
gas M gas (t) and metallicity Z g (t) as a function of time 
through our model of chemical evolution. We then use 
the scaling of the extinction curve with gas column den- 
sity and metallicity described in Guiderdoni & Rocca- 



Volmerange (1987) to compute the face-on optical depth 
of our objects at any wavelength: 



Ax 
Ay 



Zg{t) 



(JVh(*)) 



2.1 x 10 21 at em- 



ir) 



where the mean H column density (accounting for the 
presence of helium) reads: 



MgaJf) 

\A\xmpKrl 



(8) 



The second step is more delicate because it involves 
choosing a "realistic" geometry distribution for the rela- 
tive distribution of stars and dust. We model galaxies as 
oblate ellipsoids where dust and stars are homogeneously 
mixed, and scattering is taken into account. As explained 
in Paper I, the model gives a decent fit of the sample of 
local spirals analysed by Andreani & Franceschini ( 199G| ) . 

Finally, the third step involves an explicit modelling 
of the dust grain properties and sizes. We u se th e three- 
component model described in Desert et al. ( 1990| ) for the 



U / F300W, s +4 
20 25 30 




Fig. 1. Influence of the different cosmologies on the 
UV/near-IR faint counts. Dots stand for ACDM, dashes 
for OCDM, and solid lines for SCDM. Data are from Hogg 
et al. (1997) (U band), Williams et al. (1996) (F300Wab, 
B & I bands), Arnouts et al. (1997) (B band), Bertin k, 
Dennefeld (1997) (B band), Gardner et al. (1996) (B, I 
& K bands), Metcalfe et al. (1995) (B band), Weir et al. 
(1995) (B band), Smail et al. (1995) (I band), Le Fevre 
et al. (1995) (I band), Moustakas et al. (1997) (K band), 
and Djogorvski et al. (1995) (K band). 

Milky Way with polycyclic aromatic hydrocarbons, very 
small grains and big grains, and we allow a fraction of 
the big grain population to be in thermal equilibrium at 
a warmer temperature if our galaxies undergo a massive 
star burst. The weights of these four components are fixed 
in order to reproduce the relations of IR/submm colours 
with bolometric IR luminosity Ltr that are observed lo- 
cally, as detailed in Paper I. Once the full (UV/submm) 
spectral energy distributions of individual objects are 
computed following such a method, we build populations 
of galaxies for which we derive galaxy counts and redshift 
distributions. We present these results in the following sec- 
tions. 

3. The free parameters 

In addition to the cosmological parameters h, tto, Ao, £Ib 
and o-g/j-i , and to the choice of the IMF, which is assumed 
to be constant throughout a Hubble time, we basically 
have three astrophysical free parameters in the current 
version of our simple semi-analytic model: the star for- 
mation efficiency /3 , the SN heating efficiency csn, and 
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Fig. 2. Influence of the different cosmologies on the mid- 
IR/FIR/submm counts. Data are from Elbaz et al. (1999) 
(15 fim), Kawaraet al. (1998) and Puget et al. (1999) (175 
/im), Smail et al. (1997), Eales et al. (1999) and Barger et 
al. (1999a) (850/zm). Coding for the lines is the same as 
in Fig |. 



the disk truncation parameter f c which is used to compute 
the gas column density and the face-on optical depth. As 
a matter of fact, there is not much freedom in the choice 
of these parameters. 

First, the value of the star f ormation efficienc y /3 _1 de- 
duced from Kennicutt's data ( Kcnnicutt, 1998| ) is about 
(3 ~ 50 for our definition, and is valid for galaxies ranging 
from quiescent objects to very active starbursts. We refer 
to GHBM for a discussion of how this prescription actually 
compares to the data, especially to the so-called "Roberts 
times" in nearby disks, and just mention that the differ- 
ence between the value in this paper and the one used in 
GHBM ((3 = 100) stems from the different prescriptions 
used to compute disk sizes, which result in our disks being 
about twice as large as theirs. As mentioned by Kcnnicutt, 
there is a lot of scatter in the data (± 30-50 %), which, 
along with plausible systematics in the calibration of the 
different star formation estimators, should make the value 
of (3 uncertain by at least 20 %. Increasing (3 decreases 
the normalisation and slope of the optical and IR counts, 
because star formation is lower, and takes place at lower 
redshifts. 



Second, recent numerical simulations ( Thornton et al 



1998|) sjuggest that the SN heating efficiency e SN is ~ 0.09. 
However, there is much uncertainty on the actual effi- 



ciency of SN explosions in a disk galaxy, because SN bub- 
bles can blow their energy out of the disk without altering 
the cold gas (see e.g. De Young & Heckman ( |1994| ), and 
Lobo & Guiderdoni ( 1999| ) for an examination of the is- 
sue within a SAM) . Consequently, esN could be very low. 
We adopt esN = 0.03 in the following. Increasing esN de- 
creases the normalisation and slope of the optical and IR 
counts, since star formation is quenched in galaxies with 
still higher masses, that form at still lower redshifts. 

Third, the average value of the disk truncation pa- 
rameter / c , which measures the gaseous disk extension, is 
around 6, from the sample of spiral galaxies with various 
morphological types observed by Bosma ( |1981| ), and used 
by GHMB. However, this number is probably uncertain by 
about a factor 2, and it can be adjusted within this range 
in order to match the UV/optical/near-IR counts as well 
as IRAS counts, as far as this parameter fixes the amount 
of dust absorption in a disk galaxy. Increasing f c increases 
the normalisation of the optical counts and decreases that 
of the IR counts, since extinction decreases. 

As explained in Paper I, the set of Stardust spectra 
depends on (i) the mass of baryons in the galaxy, (ii) the 
star formation timescale i*, (iii) the age t of the stellar 
population, and (iv) the parameter called Jh that links 
the gas mass fraction to the gas surface density, and is 
used in the computation of the face-on optical depth. In 
the SAMs, these quantities are computed directly from the 
cold gas mass Mi, ar , the dynamical time td yn , the collapse 
rcdshift z co u and observed redshift z, and the disk expo- 
nential length rd, provided the values of the cosmological 
parameters are chosen, and the astrophysical parameters 
(3, csn and f c are fixed. Thus the implementation of the 
STARDUST spectral energy distributions within our SAM 
is very straightforward and does not bring new free pa- 
rameters. 

The above-mentioned values of the free parameters de- 
fine the so-called "quiescent mode" of star formation, sim- 
ilar to what is observed in local disks. As a reference point, 
we list, for the ACDM cosmology, the properties of a disk 
galaxy hosted by a halo of about 7.5 x 10 11 Mq, which 
collapses at a redshift ~ 3 (meaning that the age of the 
"Milky Way" -class spiral galaxy that sits in this halo is 
about 11.4 Gyr) with a spin parameter ~ 0.08. At red- 
shift 0, such a galaxy has turned about 88 % of its total 
7.5 xl0 10 M Q of cold gas into stars. Its disk exponential 
scale length is about 3.5 kpc, yielding a gaseous disk ex- 
tending to 21 kpc with an average hydrogen column of 
8 x 10 20 at cm~ 2 and a metallicity of 0.02. This, in turn 
implies a face on optical depth in the B band of 0.7 re- 
sulting in values of Mb = —19.8 and Lir — 2 x 10 10 L© 
for the face-on absolute B magnitude and the bolometric 
IR luminosity (between 3 and 1000 /im) respectively. 

Although this discussion is only valid, strictly speak- 
ing, for a given cosmology, it is unlikely that different 
cosmological parameters will significantly affect physical 
parameters like star formation or feedback efficiency. We 
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therefore consider our astrophysical parameters as inde- 
pendent of the cosmological model. In the next section, 
we keep the same values of the astrophysical parameters, 
and we study the predictions of the SAM for the sets of 
cosmological parameters that are displayed in table [i] 
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Fig. 3. Influence of the different cosmologies on multi- 
wavelength redshift distributions of galaxies. Coding for 
the lines is the same as in Fig [jj Data are from Crampton 
et al. (1995) (Canada-France Redshift Survey), and Ashby 
et al. (1996) (North Ecliptic Pole region). The predicted 
curves in the / band and at 60 (im have been renormalised 
to the total number of galaxies in the data. 



4. Sensitivity of faint counts to cosmological 
parameters 

We emphasise that our purpose here is not to determine 
values of h, f^o or Ao, but rather to answer the following 
question: what is the net effect of changing the cosmo- 
logical parameters on faint galaxy counts from the UV to 
the submm. We explore the set of cosmological parame- 
ters displayed in table [|. All the astrophysical parameters 
are fixed at the "natural values" of the "quiescent" mode 
of star formation based on the local universe. 

On Figs. |l| and |[ we show multi-wavelength counts 
obtained for the SCDM, ACDM, and OCDM cosmologies 
defined in table H. From these figures, one clearly sees 
that, in agr eemen t with Heyl et al. ( 1995| ) and Somerville 
&Primack (|2000| ), the UV to NIR counts are relatively in- 
sensitive to changes in the cosmological parameters. This 



Cosmological Mod- 
els 




Ao 


ft 






SCDM 


1.0 


0.0 


0.5 


0.015/1" 2 


0.58 


ACDM 


0.3 


0.7 


0.7 


0.015/i -2 


1.0 


OCDM 


0.3 


0.0 


0.7 


0.015 ft" 2 


1.0 



Table 1. Parameters of the different cosmologies. 



relatively "stable" behaviour extends to the FIR range. 
In sharp contrast, the differences between the predictions 
of the various cosmological models are spectacular in the 
submm range : at 850 /im and for the flux level of 10 mJy, 
the model predicts ~ 100 times more sources (or ~ 10 
times brighter sources) in the OCDM cosmology than in 
the SCDM. An interesting trend also comes out of these 
figures: with the quiescent model, any low matter-density 
universe does a better job at matching the ISOPHOT 
counts at 175 /<m and the SCUBA counts at 850 ^m than 
a critical one. Our OCDM is even able to fit the submm 
counts "naturally", without any additional ingredient. 

The influence of cosmology on the faint counts is pro- 
duced by the complicated combination of several effects. 
For instance, for lower values of the density parameter 
r2o> either with zero cosmological constant, or with zero 
curvature, we have the following changes: 

1. For halos of a given mass, the collapse occurs earlier 
on an average (see Fig A.l); 

2. The halo number density is lower (see Fig |A.l| ); 

3. With a fixed value of fls, the baryon fraction is higher, 
so that there is, on an average, more fuel for star for- 
mation per halo of a given mass; 

4. Volume elements are larger; 

5. Luminosity distances are larger; 

6. The time versus redshift relationship changes, and the 
amount of evolution undergone by the sources at any 
redshift with respect to z — is larger. 

These six effects act in different ways. If they are taken 
separately, points 3, 4, and 6 increase the slope of the 
counts whereas points 1, 2, and 5 decrease the slope. In 
addition to this, there is the effect of the fc-correction. In 
the optical, and NIR, the ^-correction is positive, and can- 
cels out partly what is occurring at high redshift in such 
a way that the faint counts are weakly sensitive to cos- 
mology. At optical wavelengths, the net effect is that faint 
counts in low matter-density universes are below those 
in the SCDM. This conclusion is opposite to the predic- 
tions of phcnomenological models based on backward evo- 
lution of the local luminosity function, under the assump- 
tion of monolithic collapse and pure luminosity evolution 



(see e.g. (Guiderdoni and Rocca-Volmerange, 1990)). The 



latter models predict that the OCDM is over the SCDM. 
The origin of this discrepancy is that the phenomenologi- 
cal models do not take points 1, 2 and 3 into account. 
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This weak sensitivity extends to the mid-IR and FIR, 
but the behaviour of the submm counts is very different. 
In contrast with the optical, NIR, mid-IR and FIR, the 
fc-correction is negative at submm wavelengths (see Fig 17 
of Paper I for an illustration), and enhances the effects of 
cosmology and evolution at high redshift. The dominant 



effects are the earlier collapse (see Fig A.l) and larger 
volume available at high redshift in low matter-density 
universes. 

Such an effect is particularly marked on the redshift 
distribution of the sources given in fig |^. These predic- 
tions arc compared with the faintest redshift survey in 
the I band (the CFRS, Lilly et al. fl!995| ), Crampton 
et al. ( ]1995| )) , with I AB < 22.5, and the North Ecliptic 
Pole Region (NEPR) survey from IRAS at 60 fim, with 
60 < S 60 < 150 mJy QAshby et al.7T99^). The CFRS sur- 



vey is correctly reproduced by the quiescent model, what- 
ever the cosmology, though all quiescent models seem to 
overpredict low-luminosity objects and produce a peak 
at too low a redshift with respect to the data (~ 0.3 in- 
stead of ~ 0.5). This is clearly due to the overproduc- 
tion of low-luminosity objects in the luminosity function. 
There is not much sensitivity to cosmology. At 60 /zm, 
the various cosmologies predict different redshift distribu- 
tions. The SDCM peaks at low redshift, without any high- 
redshift tail, as anticipated by GHBM. The ACDM and 
the OCDM peak at higher redshift, with broader distri- 
butions. The OCDM already seems to overpredict high-z 
galaxies in the NEPR at 60 /im. Finally, the sensitivity 
of the redshift distributions to cosmology is spectacular 
in the submm range. The wavelengths 175 /im and 850 
^m and the flux cuts at 100 and 2 mJy respectively cor- 
respond to on-going redshift surveys of the ISOPHOT 
and SCUBA sources. There are no firm results for these 
surveys because of identification problems (Downes et al. 
( |1999| ), Smail et al. ( 1999| )), and we prefer not to plot data. 

However, our quiescent models do not contain merg- 
ers and therefore do not account for the massive ULIRGs 
seen by ISOPHOT and SCUBA. This has to be taken 
into account in the model, though a low matter-density 
universe lessens noticeably the importance of the contri- 
bution of ULIRGs to the cosmic FIR luminosity. We now 
try to assess the impact of such mechanisms on our results 
phenomenologically. 

5. Sensitivity of faint counts to the star formation 
history 

Our simple SAM is not able to compute cither the merg- 
ing history of halos, or of the galaxies they host. However, 
we know that locally there is a tight correlation between 
major mergers on one side, LIRGs and ULIRGs on the 
other: at least 95 % of them are currently undergoing ma- 
jor mergers (see for instance Sanders & Mirabel ( 1996j )). 
It also seems fairly safe to assume that ISOPHOT and 
SCUBA sources are the high-redshift counterparts of such 



mergers. As a matter of fact, one could sum up the qual- 
itative information from currently available datasets as 
follows. First, the objects seen by SCUBA have to be 
either very massive, or very efficient to extract energy 
from the gas, simply because their bolometric luminos- 
ity is larger than 10 12 L©. Second, they have to be highly 
extinguished because most of this luminosity is emitted in 
the IR/submm. Third, for such numerous bright sources 
not to have been detected in the IRAS NEPR redshift 
survey at 60 /im, they have to be located in majority at 
redshifts greater than about ~ 1.5, which seems to be the 



case for some of the SCUBA sources flBarger et al., 1999b|) 




-5-4-3-2-10 1 2 -5-4-3-2-10 1 2 

log 10 S (Jy) log.. S (Jy) 



Fig. 4. IR counts for the SCDM "quiescent" model (solid 
lines), the "burst" model (dots), and a still more efficient 
model to extract luminosity from the gas (dashes). The 
luminosities (resp. number densities) of ULIRGs are mul- 
tiplied (resp. divided) by a factor 2 in the efficient model 
as compared to the "burst" model. 

In light of these observational facts, and as in GBHM, 
we define an ad-hoc "starburst" model, simply by push- 
ing the limits of our quiescent models (SCDM, OCDM, 
or ACDM), still powering the sources with star formation 
according to a Salpeter IMF. This consists merely in trans- 
forming a fraction of high-redshift quiescent objects into 
ULIRGs, while keeping all the parameters of the model 
fixed. The obvious interest of such an exercise is to as- 
sess whether one is able to reproduce the SCUBA source 
counts, along with preserving the quality of the fits of the 
optical counts used to calibrate the quiescent model, in 
the various cosmologies. We hereafter focus on the SCDM 
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cosmology, for which the "quiescent" mode of star forma- 
tion is unable to reproduce the submm counts. 
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Fig. 5. Multi-wavelength redshift distributions for the 
SCDM "quiescent" model (solid lines), the "burst" model 
(dots), and a still more efficient model to extract luminos- 
ity from the gas (dashes), as in Fig |[ 

In order to build such an ad-hoc model, we use the 
reasonable recipe that follows: 



1. 



A fraction of objects with halo masses larger than 
1O 12 M0 goes through a ULIRG phase when their host 
halos collapse; their SFRs and optical depths are typi- 
cally two orders of magnitude higher than those of the 
z = Milky Way (e.g. Rigopoulou et al. ( |1996[) ). We 
tune our f3 and f c parameters to obtain such proper- 
ties for the heavily-extinguished burst mode of star 
formation. Typically, we take (3 — 0.5 and f c = 0.5 for 
the starbursts. 

These ULIRGs are mainly located at z > 1.5, which 
we enforce by requiring that their fraction evolves pro- 
portionally to the squared density, i.e. as (1 + z co u) e . 
Their number density at redshift is consistent with 
the IRAS lu minosity function of ( Soifcr and Neuge- 
bauer, 199l"l). 



As a result of this phenomenological recipe, a typi- 
cal halo of mass 10 12 M Q , with reduced spin parameter 
A ~ 0.04, that collapses at redshift ~ 3, hosts by red- 
shift ~ 2.8 (180 Myr after the starburst was triggered) 
a ULIRG of size 1 kpc that has consumed 98 % of its 
10 11 Mq of cold gas initially present. The star formation 



rate averaged over this period is ~ 540 M© yr _1 . The 
starburst galaxy has a typical column density of about 
7 x 10 22 at cm~ 2 , and a metallicity of 0.03, yielding a face- 
on optical depth in the B band of 128. Its absolute B mag- 
nitude and bolometric IR luminosity (between 3 and 1000 
/im) reach Mb = —19.58 and 2 x 10 12 L Q respectively. 



U / F300W 1B +4 
20 25 




Fig. 6. 

line) . 



UV/near-IR counts for the fiducial model (solid 



Of course, such a model is quite drastic, but once 
again, it should be considered as the necessary extension 
of the quiescent models to produce the correct amount 
of FIR/submm luminosity. The interesting result is that 
such a SCDM model in which all massive objects that form 
at redshifts higher than 1.5 are ULIRGs produces almost 
enough IR/submm luminosity to match the ISOPHOT 
and SCUBA counts, as can be seen in Fig ||. This is also 
the typical luminosity one can extract from star forma- 
tion with a Salpeter IMF without ruining the UV/IR cal- 
ibration of the counts. For instance, decreasing the mass 
above which the ULIRG phenomenon occurs by an order 
of magnitude strongly decreases the optical counts. 

We also have to examine the possibility that a more 
efficient mechanism powers these sources, for instance a 
top-heavy IMF, with all the energy available through stel- 
lar nucleosynthesis being reprocessed in the IR/submm. 
The main features of such a model have been discussed 
in GHBM who take this solution to accommodate submm 
counts easily in an SCDM cosmology. We refer the reader 
to that paper for details. To test this possibility, we simply 
take our burst model and multiply the luminosity output 
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Fig. 7. Mid-IR/submm counts for the fiducial model Fig. 8. Multi-wavelength redshift distributions for the 
(solid line). fiducial model (solid line). 



of each ULIRG in the infrared per unit mass, Lir/M, 
by a factor 2, while lowering the number of ULIRGs in 
the model by 2. This is to say, we trade the number of 
sources for more luminosity per source. Fig ^ and Fig ^| 
show that the influence of such a redistribution on the 
counts is weak. Of course, any combination of luminosity 
and number density of ULIRGs is possible. 

In light of the previous work, and bearing in mind that 
we want to describe multi-wavelength galaxy counts, we 
can define a "best guess" model within a given cosmologi- 
cal model. We hereafter retain the ACDM model as a typ- 
ical example, since the optical and submm counts with the 
ACDM model and the quiescent mode of star formation 
only are intermediate between the SCDM and OCDM. We 
take h = 0.7, Q Q = 0.3, A = 0.7, Q B = 0.015fr- 2 , and 
o's.h- 1 = 1- I n terms of our astrophysical parameters, we 
keep the standard value £sn — 0.03, and we take (3 = 50, 
f c = 6 for the quiescent galaxies, and f3 = 0.5, f c — 0.5 
for the starbursts. One ULIRG dwells in each halo that 
is more massive than 10 12 M Q and collapses before red- 
shift 1.5. This ULIRG population evolves as the density 
squared at lower z, so that at redshift 0, its number density 
is about 10 -7 Mpc~ 3 , corresponding to only one ULIRG 
for 2500 halos > 10 12 Af Q . The predictions for the faint 
counts are given in Figs |^, and ^. The model provides a 
good fit of the faint counts at optical wavelengths (though 
the bright counts are slightly overestimated). The quality 
of the fit nicely compares with othe r fain t counts obtained 
from SAMs (e.g Kauffmann et al. (|l994| )). The ISOCAM 



15 /ttm data and IRAS 60 /im data are also fairly repro- 
duced, though the observed slope of the 15 /im counts 
seem to be slightly steeper than the model. The fit of the 
submm counts is also very satisfactory. 

The redshift distributions are given in Fig || The 
CFRS predictions now peak almost at the correct red- 
shift. The NEPR predictions still exhibit a high-redshift 
tail as in GHBM, in contrast with the data, but the level 
is much lower than in GHBM. We recall that the NEPR 
sample is polluted by a supercluster in the first redshift 
bin. Moreover, a recent follow-up of this sample with ISO- 
CAM at 15 /im seems to show that some of the sources 
are multiple and that the optical identifications might be 
ambiguous in these cases ( |Aussel et al., 2000 ). The rel- 
ative levels of the two peaks in the redshift distribution 
at 175 /Ltm are sensitive to the flux cut-off. Most of the 
sources in the redshift distribution for the SCUBA deep 
surveys at 850 /im are predicted to be at z > 1, but the 
comparison with data is still difficult because of identifica- 



tion uncertainties (see e.g. Barger et al. (1999b) corrected 
after Smail et al. fll999| )). 

Finally, Fig ^ shows the Cosmic Background obtained 
by integrating the faint counts, and compare the predic- 
tions with current data in the optical, IR and submm. 
Whereas introducing ULIRGs in an ad-hoc way into 
our simple models suffices to reproduce the Cosmic IR 
Background and the submm counts at 850 /mi, it falls 
marginally short of getting the required diffuse back- 
ground flux at 140 and 240 /im, though it reproduces the 



10 



J. Devriendt & B. Guiderdoni: II. Multi-Wavelength Counts from a Semi-Analytic Model 



a 



"n 1 — i 

Faint Counts 




v (GHz) 



Fig. 9. Diffuse Background light for the fiducial model 
(solid line). 



ISOPHOT counts brighter than 100 mJy at 175 /mi. These 
galaxies contribute only 10 % of the background. So this 
discrepancy may be due only to the fact that the 175 fim 
counts below 100 mJy are much steeper than our predic- 
tions. The model is too low by a factor of 2 with respect to 
the point s corrected for warm galactic dust by Lagache et 
al. ( |1999| ), which are themselves a factor of 1.5 below the 
points without such a correction by Hauser et al. ( 1998| ). 
The difficulty to fit the points might indicate that this 
correction is still underestimated. Finally, one should also 
be aware that a contribution of intergalactic dust (with 
a grey extinction curve) to the background light is also 
possible ( Aguirrc and Haiman, 1999| ). Adding these extra 
components might help reconcile models and observations. 

We conclude from these figures that this fiducial model 
gives a satisfactory estimate of the luminosity budget of 
galaxies, and allows us to interpolate or extrapolate the 
observed faint counts to other wavelengths and fainter flux 
levels. 



6. Discussion and conclusions 

In this paper, we have proposed a first implementation 
of the set of Stardust synthetic spectra into a SAM of 
galaxy formation. Although our model is quite simple, and 
cannot properly handle the merging history trees of halos 
and galaxies, we have shown that the implementation of 
Stardust is quite straightforward. This implementation 
can be easily achieved in more sophisticated SAMs. 



We have not explored a large set of statistical prop- 
erties in this paper, but have only illustrated the ability 
of this approach to reproduce the optical/IR/submm lu- 
minosity budget by producing predictions of faint galaxy 
counts at UV, visible, NIR, mid-IR, FIR, and submm 
wavelengths. As in GHBM, we have defined a quiescent 
mode of star formation which corresponds to the "natu- 
ral" values of these astrophysical parameters as they are 
suggested by local observations of disks (for [3 and f c ), or 
by numerical simulations (for csn)- 

We have then studied the influence of the cosmological 
parameters fio and A on the faint counts. In SAMs, the 
cosmological parameters influence the counts at various 
stages of the computation, through halo collapse as well 
as through the relationship of the cosmic times, luminosity 
distances and volume elements to redshifts. These quanti- 
ties intervene in the computation of the counts in different 
ways, and their effect is dimmed or enhanced by the k- 
corrections. The net result on the optical and NIR counts, 
with the influence of a positive fc-correction, is a rather 
weak sensitivity to the values of the cosmological param- 
eters. In contrast, because of the negative fc-correction 
that enhances what happens at high redshift, the submm 
counts show a strong sensitivity to cosmology. 

We know that some of the SCUBA sources detected 
at 850 /zm are the high-redshift counterparts of local 
LIRGs and ULIRGs. These objects harbour heavily- 
extinguished starbursts due to gas inflows triggered by in- 
teraction/merging. Our SAM is unable to address this pro- 
cess, and, as in GHBM, we chose to implement a heavily- 
extinguished starburst mode of star formation, by increas- 
ing the number fraction of massive objects that undergo 
this stage, as the squared density. Of course, the LIRGs 
and ULIRGs can also be powered by a top-heavy IMF. We 
have shown the strong sensitivity of the submm counts to 
the details of the ULIRG scenario, because of the negative 
fc-correction, contrasting it with the weak influence on the 
optical counts. 

We have also produced redshift distributions at vari- 
ous wavelengths and flux cuts. Here again, the sensitivity 
of the results to cosmology and evolution is dramatic at 
submm wavelengths. We got a fair fit of the CFRS red- 
shift distribution in the / band, and the redshift distri- 
bution predicted for the NEPR survey looks much bet- 
ter than in GHBM. At submm wavelengths, we predict 
a double-peaked distribution with nearby (mostly qui- 
escent) sources and distant (mostly starburst) sources. 
The relative weight of the two broad peaks is sensitive 
to cosmology and evolution. Sufficient statistics in an ob- 
servational sample could in principle help disentangling 
these effects. Unfortunately, the identification process of 
the submm sources is difficult, either in ISOPHOT or in 
SCUBA samples, and we might have to wait for multi- 
wavelength observations with forthcoming satellites such 
as SIRTF and first to solve this issue. In this context, 
our predicted counts are also a useful ingredient to anal- 
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yse current data, and to prepare observational strategies 
with these satellites. In this purview, we have proposed 
a new "fiducial model" in a ACDM cosmology which can 
take the place of the model "E" in GHBM, and which is 
available upon request. 

Merging-triggered starbursts and subsequent bulge 
formation are key-processes in the paradigm of hierar- 
chical galaxy formation. It turns out that these pro- 
cesses are particularly apparent at IR/submm wave- 
lengths, and almost invisible at optical wavelengths. So 
multi-wavelength observations are required to constrain 
the history of galaxy formation in the quiescent and star- 
burst modes. The complete merging history of galaxies 
has to be followed in detail, especially during the short 
periods of interaction and merging which produce IR lu- 
minous sources. A hybrid method of galaxy formation us- 
ing high-resolution N-body simulations to plant semi- 
analytic galaxies is being developed and will be used to 
qua ntify the importance of merging processes (Hatton et 
al. ( 2000 ) and following papers). In this context, the reader 
is invited to notice the following point: any change in the 
luminosities and number densities of the ULIRGs results 
into spectacular changes in the submm counts, because 
of their super-Euclidean regime. This sort of "instabili- 
ties" in the behaviour of the faint counts is going to be 
very useful to constrain the luminosity and duration of the 
starbursts in more refined models. 

Throughout this paper, we have assumed that dust 
heating is powered by star formation. The other possible 
engine is the presence of heavily-extinguished AGNs lo- 
cated at the centers of ULIRGs. In local samples, heating 
is dominated by starbursts, except in the most luminous 
galaxies (Genzel et al. ( |1998| ), Lutz ct al. ( |1998| )). So far, 
wc do not know how this situation evolves with redshift. 
The redshift distribution of the SCUBA sources seems to 
peak around redshift 2.5 where the obser ved redshi ft dis- 
tribution of quasar activity may also peak ( Pei, 1995 ). Per- 
haps a sophisticated combination of AGNs and starbursts 
with a top-heavy IMF is responsible for the evolution seen 
by SCUBA, and we signal that attempts to include both 
types of ingredients are already discussed in the literature 
(e.g. Blain et al. (1999)). However, properly disentangling 
all these components definitely requires a more sophis- 
ticated model than the one presented here. We therefore 
defer this issue to a later study, but remark that AGNs can 
be implemented self-consistently within the framework of 
the hybrid method ( Kauffmann and Haehnclt, 2000| ), so 
that one could hope to set quantitative limits on their 
respective contributions to the submm fluxes. 

The situation is indeed as complex from the theoret- 
ical as well as the observational points of view. That is 
why this simple modelling should be considered as an ex- 
ploratory step to probe the complex issues involved in 
multi-wavelength counts, and to design a more satisfying 
approach of the problem. 
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Appendix A: Dark matter halos in any cosmology 

We suppose that perturbations of the matter density field 
when the universe becomes matter dominated are com- 
pletely characterised by their power spectrum P(k): 



P(k) oc fc"T 2 (fc), 



(A.l) 



where T(k) is the transfer function (see fit given in Ap- 
pendix G of Bardeen et al. ( 1986| )). We further assume a 
post-inflation Harrison-Zel'dovich power spectrum (n = 
1) for these perturbations, and take the shape parameter 
T used in the computation of T(k) to be flSugiyama, 1995| ): 



r = ilo^iexp 



-n B (i + 



r» ' 



(A.2) 



where tt is the current matter density (in critical 
density units), £Ib is the baryon density, and h — 
i?o/(100 km/s/Mpc) is the reduced Hubble constant. 

In the linear regime, the equation of motion is solved 
for the expansion factor, a, and the solutions for the 
growth of the density contrast 5 (see e.g. Peebles (1980)) 
are derived. There are two such solutions (Heath, 1977), 



which form a complete set ( ^el'dovich, 1965| ) and read 
D d [z] = H {n (l + zf + (1 - O - A )(l + zf + ^A}3) 

(A.4) 



D g [z] = D d [z] 



(1 + x) 



D 3 \ 



dx, 



where the subscripts d and g respectively stand for the 
decaying and growing modes, and Ao = A/ (3Hg) is the re- 
duced cosmological constant. If one further assumes that 
the initial peculiar velocity of the perturbation is zero (i.e. 
that the perturbation simply moves along with the ex- 
panding universe), the density contrast 8[z] in the linear 
regime grows as: 



5[z] 



|fio(l 



H, 



l-Sh-^o> -rDg[z]Si, (A.5) 



where the subscript i stands for the initial quantities. 

In the non-linear regime, one considers an isolated 
spherical perturbation of radius r^, at time ti, which has 
a uniform overdensity Si = (p p [U] — pb[U]) / Pb[U] with re- 
spect to the background (Si <C 1), and encloses a mass 
M = 4/3Trrfpb[ti](l + Si). We assume that the perturba- 
tion is bound, and that its peculiar velocity is nil. The 
equation of motion is integrated to obtain the time t m at 
which the perturbation reaches its maximum expansion 
radius r m : 



t, 



U + 



,1/2 



dx, 



(A.6) 



where 



D p [x] = X x 3 + {(1- Oq- A )(l + Zi f (A.7) 
- n S l {l + z l ) 3 }x + n (l + S l )(l + z l ) 3 , (A.8) 

and r m /ri is the first real r oot > 1 of the cubic equation 
(c.f. Richstone et al. (|l992|) ): 

D p [x] = 0. (A.9) 

After it has reached this maximum radius at time t m , 
the perturbation, by symmetry, collapses on a time scale 
tcoii = 2£ m — ti. Thus, with the previous equations, the 
critical density contrast So linearly extrapolated till today 
(eq. AJ3) is explicitly related to the collapse redshift z co u 
of the perturbation for any cosmology, just by computing 
the redshift to which the collapse time of a perturbation 
with overdensity Si corresponds in the unperturbed uni- 
verse (provided t co u is smaller than the age of the uni- 
verse): 

tcoll '■■ 



1 



con 0- + x)D d [x 



dx. 



(A.10) 
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Fig. A.l. Evolution of the formation rate of dark mat- 
ter halos (masses indicated on the figure) for three differ- 
ent cosmologies: SCDM (solid line), OCDM (dashes), and 
ACDM (dots). The parameters used for these models are 
given in table |[ 

If the resulting virialized perturbation can be approxi- 
mated by a singular isothermal sphere truncated at virial 
radius r V i r , then ?w/r m is th e solu tion of the following 
cubic equation (see Devriendt (|1999| )): 



3GM y 5 5GM ' v ' 

Finally, the velocity of a test particle moving on a circular 
orbit in the isothermal sphere reads: 



r- GM 2A TT9 , 
V d = V2<t = J ^alrl 



-o'vir- ( A - 12 ) 

Implementing these results into the peaks formalism 
described in GHBM enables one to derive formation rates 



for dark matter halos as a function of redshift. Fig A.l 
illustrates this halo formation rate for three typical cos- 
mologies gathered in table |l|. 



